Constraints on Spin-Dependent Short-Range Interaction between Nucleons 
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We report on the search for a new spin-dependent P- and T- violating interaction between nucleons 
mediated by light, pseudoscalar bosons such as the axion which was invented to solve the strong CP 
problem. Our experimental approach is to use an ultra-sensitive low-field magnetometer based on 
the detection of free precession of co- located 3 He and 129 Xe nuclear spins using SQUIDs as low-noise 
magnetic flux detectors. In the presence of an unpolarized mass the precession frequency shift was 
measured to determine the coupling of pseudoscalar particles to the spin of the bound neutron. For 
the force range from 3-10 -4 m - 10 _1 m, corresponding to the preferred mass range between 2 /xeV 
and 500 fieV for the axion or axion-like particles, we improved the upper bounds on the products of 
scalar and pseudoscalar couplings (g s g P ) by up to 4 orders of magnitude. 
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Axions are light, pseudoscalar particles that arise in theo- 
ries in which the Peccei-Quinn U(l) symmetry has been 
introduced to solve the strong CP problem pQ. They 
could have been produced in early stages of the Universe 
being attractive candidates to the cold dark matter that 
could compose up to ~l/3 of the ingredients of the Uni- 
verse [2]. Several constraints from astrophysics, cosmol- 
ogy, and laboratory experiments have been applied in 
order to prove or rule out the existence of the axion, i.e., 
constrain the axions mass m a and/or its couplings. The 
mass range, in which axions are still likely to exist, could 
thus be narrowed down to a window reaching from fieV 
[3] up to some meV [4] (axion window). 
Most axion searches look for the conversion of a galactic 
[5], solar [6], or laboratory [7] origin axion into a pho- 
ton in the presence of a static magnetic field. However, 
any axion or axion-like particle that couples with both 
scalar and pseudoscalar vertices to fundamental fermions 
would also mediate a parity and time-reversal symmetry- 
violating force between a fermion / and the spin of an- 
other fermion f a , which is parameterized by a Yukawa- 
type potential with range A and a monopole-dipole cou- 
pling given by [8]: 

a is the spin vector and A is the range of the Yukawa-force 
with A= h/(m a c). Thus, the entire axion window can 
be probed by searching for spin-dependent short-range 
forces in the range between 20 /im and 0.2 m. gf and 
gf* are dimensionless scalar and pseudoscalar coupling 
constants which in our case correspond to the scalar cou- 
pling of an axion-like particle to a nucleon (g£ = gjf) 



and its pseudoscalar coupling to a polarized bound neu- 
tron (gf* = g™). r is the unit distance vector from the 
bound neutron to the nucleon. The potential given by 
Eq. [I] effectively acts near the surface of a massive unpo- 
larized sample as a pseudomagnetic field and gives rise to 
a shift A^ sp = 2-Vs/ft, e.g., in the precession frequency of 
nuclear spin-polarized gases ( 3 He and 129 Xe), which ac- 
cording to the Schmidt model [9 can be regarded as an 
effective probe of spin-polarized bound neutrons. The po- 
tential Vs is obtained by integration of V sp (r) from Eq.[l] 
over the volume of the massive unpolarized sample aver- 
aged over the volume of the polarized spin-sample. Based 
on the analytical derivation of Ve j00 f° r disc-shaped spin- 
and matter samples with respective thicknesses D and d 
[TP] , we can derive the following expression for V% given 

by 

= Fs.oo • >7(A) = 2ttNk^ ■ e~ Ax ^ x 

(l _ e -J>/A) . (i _ e -<VA) . 77(A) . (2) 

77(A) takes account for the finite size in transverse di- 
rection of our cylindrical samples and Ax represents the 
finite gap between them. Furthermore, k = h 2 'g^ g™ / '(87r • 
m n ) and TV is the nucleon number density of the unpo- 
larized matter sample. >7(A^ is determined numerically 
for our cylindrically shaped spin- and matter samples at 
"close" -position (see Fig.[T]). 

Our experimental approach to search for non-magnetic, 
spin-dependent interactions is to use an ultra-sensitive 
low-field comagnetometer based on detection of free spin 
precession of gaseous, nuclear polarized samples [TT] . 
The Larmor frequencies of helium and xenon in a guid- 
ing magnetic field B are given by o; L ,He(Xe) = 7He(Xe) ■ B, 
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1 77(A) can be expressed reasonably well by the fit function 77^ (A) = 
(1 + 27.8- A 1 - 31 )/(l + 234 -A 1 - 31 ). 
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FIG. 1. Sketch of experimental setup. The dewar houses a 
total of 304 low-T c DC-SQUID magnetometers which are di- 
vided up into 19 identical modules. The lower plane SQUIDs 
in module D, E, and I marked as (•) are used to detect the 
3 He/ 129 Xe free spin precession. The center of the cylindrical 
spin sample cell (D — 60 mm, 0d = 58 mm) has an average 
distance of z = 66 mm to the sensors. The cells relative posi- 
tion in (x,y)-plane projection is indicated by the grey square. 
SQUIDs in module S marked as (■) are used for the gradio- 
metric sensor arrangements. The unpolarized mass (cylindri- 
cal BGO crystal: d = 70 mm, (dd — 60 mm) can be moved 
along the x-axis (5-field axis) to "close" (Ax c = 2.2mm) and 
"distant" (Axd = 170mm) position and vice versa (see text). 
This is accomplished by a piston driven glass tube with the 
BGO fixed at its cell-facing side. The two arrangements left 
(C) and right (R) are used. 

with 7He(Xe) being the gyromagnetic ratios of the respec- 
tive gas species [H [13] with 7 He /7x e = 2.75408159(20). 
The influence of the ambient magnetic field and its tem- 
poral fluctuations cancels in the difference of measured 
Larmor frequencies of the co-located spin samples 

AuJ = CJ H e - — * ^Xe • (3) 

7Xe 

On a closer look, a resulting constant frequency shift 
AcJii n , e.g., due to Earth's rotation, is not compensated 
by co-magnetometry. That is discussed in [14 , together 
with frequency shifts due to the generalized Ramsey- 
Bloch-Siegert shift which are directly proportional to the 
particular net magnetization A He ( Xe ) ■ e -t / T2 > He ( Xe ) and 
which will show up in the weighted frequency difference 
Acj(t), too: 

_ t _ t 

Au(t) = Au; lin +e H e-^He-e T2 *' He +e Xe -^Xe-e T2 * Xe • (4) 

Accordingly, its equivalent, the weighted phase difference 
A$(t) = $He(£) - ^'$xeW, is sensitive to a phase drift 
given by 

_ t 

A$(t) = $ + Au; lin • t - e He • T 2 * He • A Ue • e T ^ 

-exe-T* Xe -A X e-e"^ . (5) 

Due to the knowledge of these side effects, any anoma- 
lous frequency shifts generated by non-magnetic spin in- 
teractions, such as the quested short range interaction, 



can be analyzed by looking at Au(t) and A3>(£), respec- 
tively. A sudden frequency change Au; sp stemming from 
the pseudoscalar Yukawa potential V sp (r) would occur 
at an instant t = to, e.g., by moving a massive matter 
sample close to the precessing nuclei. This would lead 
to an additional linear phase drift Au sp • t in Eq.[5]for 
t > to. For further analysis, it is useful to develop Eq.[5] 
in a polynomial of 5th ordei[^] around to. The weighted 
phase difference A3>(£) can then be described by 

A^ fit (t f ) = a + b(t') -t' + c-t^ + d-t'^ + e-t'^ + f-t'^ (6) 

with t' = t — to. The coefficient of the linear term now 
reads 

b(t') = Auj hn + Auj sp (t') - e He • A Ue - e Xe ■ A Xe . (7) 

Note that Auj sp (t') = 2tt • A^ • 6(±£ / ) F]is the only time 
dependent term in Eq.[7| so that a cnange of b(t') at 
t = to would directly indicate the existence of the short 
range interaction. This statement, the fit-function of 
Eq.[6]is based on, is valid as long as susceptibility related 
artefacts of the BGO crystal at its respective positions 
do not lead to different relaxation times T|. With our 
special choice of t' = t — to, the linear coefficient of the 
Taylor expansion does not dependent on T 2 * and thus 
is insensitive to possible changes in T|. The impact 
of the T 2 * -dependence of higher order terms on b(t') is 
discussed in detail in section systematic uncertainties. 

The experiments were performed inside the magnetically 
shielded room BMSR-2 at the Physikalisch-Technische 
Bundesanstalt Berlin (PTB)[15j. A homogeneous guid- 
ing magnetic field of about 350 nT was provided inside 
the shielded room by means of a square coil pair (B x - 
coils) of edge length 1800 mm. A second square coil pair 
(By-coils) arranged perpendicular to the I^-coils was 
used to manipulate the sample spins, e.g., tt/2 spin flip by 
non-adiabatic switching [TT] . The main parts of the ex- 
perimental setup within BMSR-2 are shown in Fig.[l] For 
detection of the spin precession we used a multi-channel 
low-T c DC-SQUID device which was originally designed 
for biomagnetic applications [16j[I7]. The 3 He/ 129 Xe nu- 
clear spins were polarized outside the shielding by means 
of optical pumping. Low-relaxation cylindrical glass cells 
(GE180) were filled with the polarized gases and placed 
directly beneath the dewar as close as possible to the 
SQUID sensors. The SQUID sensors detect a sinusoidal 
change in magnetic flux due to the spin precession of the 
gas atoms in the glass cell. In order to obtain a high com- 
mon mode rejection ratio, three first order gradiomet- 
ric sensor combinations were used in order to suppress 



2 The criterium to choose a polynomial expansion up to the 5th 
order was that the reduced x 2 /d.o.f. of the fit equals 1. 

3 (±) in the argument of the Heaviside step function has to be set 
(-) for the sequence c— ^d and (+) for the reverse one d— >-c. Fur- 
thermore, for runs j = 1, 2, 3 the BGO was moved at fo=8700 s, 
otherwise at to =10800 s . 
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environmental disturbance fields like vibrational modes. 
Fig. [I] shows their positions with respect to each other 
and with respect to the 3 He/ 129 Xe sample cell. The sys- 
tem noise of the SQUID gradiometer configurations was 
between 3 fT/v / Hz and IOiT/a/Hz in the range of the 
3 He/ 129 Xe spin-precession frequencies, i.e., 4 Hz < vl < 
12 Hz, while typical signal amplitudes reached 10 pT for 
helium and 3 pT for xenon at the beginning of the spin 
precession cycle. The optimum conditions in terms of 
long transverse relaxation times (T 2 *) and high signal-to- 
noise ratio (SNR) were met at a gas mixture with pres- 
sures of 3 He : 129 Xe : N 2 « (2 : 8 : 35) mbar, typically. 
Nitrogen was added to suppress spin-rotation coupling 
in bound Xe-Xe van der Waals molecules [181 HI]- As 
unpolarized matter sample we used a cylindrical BGO 
crystal (Bi4Ge30i2, p=7.13 g/cm 3 ). BGO has a high nu- 
cleon number density, is a non-conductive material that 
shows low Johnson-Nyquist noise and is said to have an 
unusual magnetism-related behaviour in weak constant 
magnetic fields (Xmag ~ ppm) [20 - 22 . For systematic 
checks, primarely to look for magnetic susceptibility re- 
lated false effects, the BGO crystal could be placed left 
(£) and right (3?) with respect to the 3 He/ 129 Xe sample 
cell (see Fig. IT]). Since V sp (r) oc a • f, A^ sp changes its 
sign in going from C to That has to be considered 
by averaging the C and 5ft results. On the other hand, 
Az/ sp drops out averaging C and 5ft without sign change. 
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FIG. 2. (a) Extracted frequency shifts Ai/ sp (with correlated 
la error) of the 10 measurement runs. The triangles spec- 
ify the 5ft, the circles the C arrangement of the BGO crystal. 
Full symbols indicate the c— »d sequence, hollow symbols the 
opposite case (d— >-c). (b) Results Azy c h ec k obtained from the 
LV-data using the same fit-model (Eq.[6]). Since no mass was 
moved, we expect no shift in the spin precession frequency. 
The rightmost symbols in both plots (stars) indicate the re- 
spective weighted means (la error). 



is then extracted from 



Av sr> — — 

sp 2tt • (1 



b d 



7He \ 

7Xe ' 



(9) 



For 6 runs (# : 2, 3, 4, 5, 9, 10), the BGO crystal 
was positioned at C otherwise at 5ft. For all C runs, 
In case of a non-zero spin-dependent axion fermion inter- tne results are multiplied by (-1): A^ sp = -A^ sp? £. In 



action., a shift Az/™ in the weighted frequency difference 



(Eq.[3| can be extracted from respective frequency mea- 
surements in "close" and "distant" position given by 



A^ 



2VS 



7He 
7Xe 



(8) 



assuming Vs,H e = ^s,Xe = Vs,n = Ve (Schmidt model) 



and 



We performed 10 measurement runs lasting approxi- 
mately 9 h each. For each measurement run after to « 3h 
the BGO crystal was moved from "close to distant" posi- 
tion (c— >-d) or vice versa (d— >*c). The asymmetric timing 
takes account for the smaller SNR in the second measure- 
ment block due to the exponential damping (T 2 * ) of the 
signal amplitude which was T 2 * He « 48 h and T 2 * Xe « 5h, 
typically. By this measure, comparable statistics was ob- 
tained for both BGO positions. 

As discussed in detail in [13] . the data of each run were 
divided into sequential time intervals of r = 3.2 s. For 
each obtained sub-data set, a x 2 - mn iimization was per- 
formed using an appropriate fit-function to extract the 
phases (/>H e , </>x e and the frequencies uo^ ei uJxe with the 
corresponding errors. In a further step, the accumulated 
phase $He(Xe)(^) was determined for each run in order 
to derive the weighted phase difference A<£ (t). Then we 
fitted Eq. [6] simultaneously to the data of A$(t) deter- 
mined for the three gradiometers for each measurement 
run. From the resulting fit parameters a, 5 C , b<i, c, e, 
/ and by use of Eqs. [8]and[7| the frequency shift Az/ sp 



Fig. values A^ sp for the individual runs are shown 
together with their correlated la errorQ Calculating 

the weighted mean we get Av sp = (—2.9 ± 2.3) nHz. As 
consistency check, we re-analysed our 2009 data, where 
we looked for a possible Lorentz-violating (LV) sidereal 
frequency modulation [13] . Since no mass was moved, 
b c = bd should hold, using the fit-function of Eq.[6]and 
a hypothetical time to = 10800 s. Fig. shows the re- 
sults A^ c heck for all 7 measurement runs, again with their 
correlated la errors. The weighted mean of the LV-data 
gives A^check = (-1.4 ± 3.4) nHz. The x 2 / d -°- f of the 
data to their respective weighted means (Ai/ sp , Az/ C h ec k) 
gives 2.29 and 2.38, indicating that the errors on the 
measured frequency shifts (Fig. ^) are somewhat under- 
estimated. In order to take this into account, the errors 
were scaled to obtain a x 2 /d.o.f of one, as recommended, 
e.g., by [H [24]. 

At the 95% C.L., our results for the measured frequency 
shifts are: 



A^ sp = (-2.9 ± 6.9) nHz 
A^check = (-1.4 ±10.5) nHz 



(10) 
(11) 



4 The correlated errors are calculated as square root of the diagonal 
elements of the covariance matrix of the least x 2 fit-model of 
Eq.|6]with the proper statistical weights. The uncorrelated errors 
are about a factor of 30 smaller and not included in the error bars 
shown in Fig. ^] 
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indicating that i) we find no evidence for a pseudoscalar 
short-range interaction mediated by axion-like particles 
and ii) the cross check analysis of our LV-data within 
the error bars is also compatible with zero, as expected. 

Discussion of systematic uncertainties: 
The movement of the BGO-crystal can produce corre- 
lated effects that may mimic a pseudoscalar frequency 
shift or even compensate the effect we are looking for. 
Two effects caused by a non-zero magnetic susceptibilty 
of the BGO have to be considered. 

a. ) The BGO at "close" position slightly changes the 
magnetic field across the volume of the 3 He/ 129 Xe sam- 
ple cell. This effect drops out to first order due to comag- 
netometry. To second order, however, the imperfect geo- 
metrical overlap of helium and xenon resulting from the 
difference in their molar masses leads to a difference (Az) 
in their center of masses (c.o.m) in the gravitational field 
of the earth according to the barometric formula. For 
our cylindrical sample cell this value is Az = 1.2- 10 _7 m. 
With Xbgo = — 19 ppm, which is the susceptibility of the 
BGO in the high field limit (B > 0.1 T) [2T], a conserva- 
tive estimation of the systematic effect can be calculated. 
Therefore the field distribution across the 3 He/ 129 Xe 
sample cell, with BGO in "close" position and B x = 
350 nT, was calculated with COMSOL Multiphysics, a 
finite element analysis software. From that we derived 
induced magnetic field gradients in vertical direction of 
\(dB/dz)\ ind < 0.08pT/cm, which lead to a frequency 
shift of Az/™ = Az-\(dB/dz)\ ind - 1U e/27r = 0.03 nHz. 
Compared to the measured frequency shift (Eq.[l0|), this 
systematic effect is negligible. 

b. ) More serious is the fact that a change of the magnetic 
field gradient by the BGO also influences the T 2 * -times of 
3 He and 129 Xe. According to [HJ [25J [26] , the change in 
T 2 * for a spherical cell of radius R can be expressed by 



AT 2 * 



(T 2 *) c - (T 2 *) d 
(T 2 *)2 



8i? 4 7 2 
175D GM 



({{dB x /dx) F + (dB x /dx) ind ) 2 c - ((dB x /dx) F ) 2 d ) , 

(12) 

where D GM is the diffusion coefficient of the gas mixture 
[IB with L>g e M « 17cm 2 /s and L>£ e M « 4.7cm 2 /s, respec- 
tively. For our cylindrically shaped cell we use here the 
approximation D /2 = D/2 w R = 30 mm. The BGO 
as diamagnet shows its strongest influence ("close" po- 
sition) along the direction of the applied B x -fi.e\d which 
was calculated to be (dB x /dx) ind = ±0.47pT/cm across 
the sample cell. The (+) and (-) signs indicate the two 
BGO arrangements C and 5ft, respectively. 
Knowing (dB x /dx) F of the guiding field inside BMSR- 
2, AT 2 */(T 2 *) 2 can be calculated, in principle. The more 
direct approach, however, is to extract T 2 * and thus AT 2 * 
via the exponential decay of the signal amplitudes with 
the BGO in "close" and "distant" position. The most ac- 
curate distinction between (T 2 *) c and (T 2 * ) d was obtained 



through a fit to the amplitude ratio Ax e {t f )/^He(t f ) 
given by 



/fit (tO = W ■ e"*'/^ 



withT e * ff = T* He .T* Xe /(T* He 



T* Xe ). Using 



AT e * ff 

(t: s ) 2 



AT 2 * He 



AT 2 * Xe 



-0.5- 



ATr 



2, He 



(5,He) (T$,xe) ( T lu e y 



(13) 



(14) 



(see Eq.[l2]), a possible T 2 *-change of | AT 2 * He | < 160s can 
be inferred from this fitting procedure. From that the 



systematic frequency shift A^ sy 2 s on b(t f ) due to higher 
order terms of Eq.[6j which has the same signature as 
the pseudoscalar phase shift, can be estimated to b^] 



^sys 



< 



(^He 



(7Ixe) £ 



2tt(1 - 7He/7Xe) 



(15) 

With ^HetXe) = e H e(Xe) ' ^He(Xe) ' T 2 *He(Xe)' Here We USed 

Eq.[9j replacing b c and bd by the temporal means 2 • c c • 
{t') t , and 2-Cd' (t') t , of the quadratic term (Eq.|6| with 
c c(d) = {(£ He /2)/(T* He ) 2 + (^xe/2)/(T 2 * Xe ) 2 } c(d) . 
From the fit function (Eq. [5| applied to the data, values 
for the respective £?He(Xe) amplitudes can be extracted, 
i.e., -25 < En e < -4 and Ex e ~ 0.2. Finally, taking 
(t') t , ~ to/2 and by use of Eq. 14 we get 



Aw 



sys 



AT* 



(^,He) 2 



1^2* He + 2T« X J ^ 



< v 2 ' He/ ; x — , w 0.2 nHz 



2tt(1 - 7He/7Xe) 



(16) 

From that, a conservative estimate on the systematic er- 
tot can be made with |A^ sy 2 s | 
That brings us to the final result 



±0.4 nHz (95% C.L.). 



A^ sp = (-2.9 ± 6.9 ± 0.4) nHz (95% C.L.) 



(17) 



for the measured pseudoscalar frequency shift. 

From the total error S(Au sp ) = ±7.3 nHz we can then 

derive exclusion bounds for |<7^#p| using Eq.|2] and 

|5(Ai/ sp )| >2-V£/h which are shown in Figji] 
We have substantially improved the bounds on a spin- 
dependent short-range interaction between polarized 
(bound) neutrons and unpolarized nucleons over most 
of the axion window, tightening existing constrains on 
axion-like particles heavier than 20 /ieV by up to four 
orders of magnitudes. 

And there are clear strategies on how to improve our 
experimental sensitivity: i) Close contact of the spin 



Contributions of terms n > 3 (Eq.[6]) to the systematic error can 
be neglected. 
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system with the matter sample. For Ax ~ mm, our 
present measurement senitivity will significantly increase 
for A < 10 -3 m (see Figj3|. ii) Moving the spin/matter 
sample more frequently between its set positions (c<H>d 
and/or C <H> 51). This results in a different time structure 
for the linear term in the fit model of Eq. [6] such that the 
correlated error approaches the uncorrelated one. This 
was demonstrated in [14], already, iii.) Magnetic suscep- 
tibility related artefacts have to be eliminated by tak- 
ing zero-susceptibility matched matter samples (Xmag ~ 
ppm) as it is common practice in high resolution NMR 
spectroscopy [27] . 



X(m) 



FIG. 3. The experimental 95% confidence upper limit on 
\9s^ 9p | plotted versus A, the range of the Yukawa-force with 
A = h/(m a c). The axion window is indicated by the light 
grey area. (1): result of [28], (2): result of [29], (3): result of 
[30], (4): result of [31], (5): result of [32], (6): result of [33], 
(7): this experiment. See [3l] for bounds on the pseudoscalar 
short-range force between polarized electrons and unpolarized 
nucleons. The expected results for Ax^ mm (8) using the 
same data set demonstrates the gain in measurement sensi- 
tivity for A < 10 -3 m by reducing the minimal gap between 
the polarized and unpolarized matter sample. 



This work was supported by the Deutsche Forschungsge- 
meinschaft (DFG) under contract number BA 3605/1-1 
and the research center "Elementary Forces and Mathe- 
matical Foundations" (EMG) of the university in Mainz, 
and by PRISMA cluster of excellence at Mainz. We are 
grateful to our glass blower R. Jera for preparing the low 
relaxation glass cells from GE180. 



[1] CD. Peccei, and H.R. Quinn, Phys. Rev. Lett., 38, 1440 
(1977) 

[2] J. Jaeckel, and A. Ringwald, Annu. Rev. Nucl. Part. Sci., 
60, 405 (2010) 

[3] S. Asztalos et al., Phys. Rev. D, 69, 011101(R)) (2004) 
[4] G. Raffelt, Annu. Rev. Nucl. Part. Sci., 49, 63 (1999) 
[5] S. Asztalos et al., Phys. Rev. Lett., 104, 041301 (2009) 
[6] E. Arik et al., J. Cosmol. Astropart. Phys., 02, 008 (2009) 
[7] K. Ehert et al., Phys. Rev. Lett., 689, 149 (2010) 
[8] J.E. Moody, and F. Wilczek, Phys. Rev. D, 30, 130 
(1984) 

[9] Th. Schmidt, Zeitschrift fiir Physik A Hadrons and Nu- 
clei, 106, 358 (1937) 
[10] O. Zimmer, Phys. Lett. B, 38, 685 (2010) 
[11] C. Gemmel et al., Eur. Phys. Journal D, 57, 303 (2010) 
[12] M. Pfeffer and O. Lutz, J. Magn. Res. A, 108, 106 (1994) 
[13] International council for Science: Committee on Data 
for Science and Technology (CODATA), www.codata.org 
(2007) 

[14] C. Gemmel et al., Phys. Rev. D, 82, 111901 (2010) 
[15] J. Bork et al., Proc. Biomag., 2000, 970 (2000) 
[16] D. Drung, Physica C, 368, 134 (2002) 
[17] M. Burghoff et al., IEEE Trans. App. Supercon., 17, 846 
(2007) 



[18] B. Chann et al., Phys. Rev. Lett., 88, 113201 (2002) 

[19] B.C. Anger et al., Phys. Rev. A, 78, 043406 (2008) 

[20] E.A. Kravchenko, V.G. Orlov, and M.P Shlykov, Russian 

Chemical Reviews, 75(1), 77-93 (2006) 
[21] S. Yamamoto, K. Kuroda and M. Senda, IEEE Transac- 
tions On Nuclear Science, 50, 1683 (2003) 
[22] B. C. Grabmaier, and R. Oberschmid, Phys. Stat. Sol., 

(a) 96, 199 (1986) 
[23] The Particle Data Group, Phys. Lett. B, 592, 14 (2004) 
[24] W.H. Press, S.A. Tenkolsky, W.T. Vetterling, and B.P 
Flannery, Numerical Recipes 3rd Edition: The Art of 
Scientific Computing., Cambridge University Press, New 
York (2007) 

[25] G.D. Cates, S.R. Schaefer, and W. Happer, Phys. Rev. 

A, 37, 2877 (1988) 
[26] W. Kilian et al., Eur. Phys. J. D, 42, 197 (2007) 
[27] R. Ravi et al., J. Magn. Res., textbf205, 63 (2010) 
[28] S. Baeler et al., Phys. Rev. D, 75, 075006 (2007) 
[29] T. Jenke et al., arXiv: 1208.3875vl (2012) 
[30] A.P Serebrov et al., JETP Letters, 91, 6 (2010) 
[31] A.N. Youdin et al.,Phys. Rev. Lett., 77, 2170 (1996) 
[32] M. Bulatowicz et al., arXiv: 1301.5224vl (2013) 
[33] P-H. Chu et al., Phys. Rev. D, 87, 011105(R) (2013) 
[34] S.A. Hoedel, Phys. Rev. Lett., 106, 041801 (2011) 



